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In  recognizing  a  grasped  object  or  grasp  an  object  stably  with  a  multifingered 
robotic  hand,  tactile  sensors  mounted  on  robot  fingers  have  been  identified  as  es¬ 
sential  detecting  tool.  In  general,  grasp  stability  has  two  requirements:  sliding 
avoidance  and  excess  force  avoidance.  Hence,  it  is  required  that  a  tactile  sensor  be 
able  to  detect  the  nature  of  the  force  distribution  which  exists  between  the  surface 
of  the  grasped  object  and  the  robot  finger.  - - - — 


Recently^much  work  has  been  devoted  to  the  tactile  sensing  problem.  Several 
kinds  ofjthe  tactile  sensors  have  been  designed  based  on  electro-optics,  piezoresistive, 
or  piezoelectric  properties, etc..  Typically,  these  sensors  are  not  used  to  measure  the 
contact  force  directly,  but  to  measure  the  interior  strain  or  stress  in  an  elastic 
finger  pad.  When  sharp  and  rigid  objects,  e.g.  a  wedge  indentor,  indent  an  elastic 
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material,  very  high  stress  is  developed  at  the  contact  surface.  This  stress  is  reduced 
by  distance  from  the  contact  area,  so  a  fragile  sensor  would  be  better  protected 
within  a  layer  soft  material.  Another  advantage  of  using  an  elastic  covering  is  that 
the  contact  area  becomes  large  so  that  grasp  stability  may  be  enhanced  and  the 
features  of  grasped  objects  better  distinguished. 


The  goal  of  this  paper  is  to  study  how  a  surface  force  profile  may  be  estimated 
from  the  information  on  strain  or  stress  distribution  detected  by  tactile  sensors. 
This  problem  is  referred  to  as  the  inverse  problem  since  we  can  consider  the  stress 
or  strain  within  an  elastic  material  as  the  response  due  to  surface  loading, 
inverse  problem  can  be  treated  as  a  deconvolution  problem  or,  more  generally, 
the  problem  of  solving  an  operator  equation  of  the  first  kind.  It  is  of  interest  to 
determine  how  such  an  ill-posed  problem  may  be  solved  using  appropriate  regular- 
ization  (Tikhonov  1977).  In  this  paper,  we  will  not  consider  the  effects  of  the  noises 
(or  suppose  the  observations  have  been  passed  through  some  kind  of  filter)  and  only 
pay  attention  to  deriving  particular  operators,  analyzing  their  properties  and  using 
the  Discrete  Fourier  Transform(DFT)  approach  to  solve  the  associated  equations. 

In  Section  2,  we  will  study  a  particular  model  of  the  elastic  material,  i.e.  the 
relation  between  the  surface  load  and  stress  or  strain  beneath  the  surface.  In  Section 
3,  we  will  investigate  an  approximate  method  to  solve  the  inverse  problem  for  some 
special  cases.  In  Section  4,  we  will  give  some  examples  and  compare  the  results. 


We  are  collaborating  closely  with  a  research  group  at  the  Naval  Research  Lab¬ 
oratory  under  Dr.  M.  Peckerer,  involved  in  the  fabrication  of  silicon-based  tactile 
sensor.  We  are  also  comparing  our  numerical  results  against  analog  network  solu-)n 
tions  to  the  inverse  problem  of  this  paper. 
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2.  The  model  of  the  elastic  material 
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In  this  section,  we  attempt  to  find  the  relation  between  the  load  profile,  which 
exists  on  the  surface  of  a  half  space  of  elastic  material,  and  the  distribution  of  the 
stress  or  the  strain  beneath  the  surface  of  that  material.  To  specify  the  load,  we 
shall  just  consider  the  contact  force  profile  rather  than  the  contact  displacement 
profile  for  the  following  reasons: 

•  The  contact  stress  is  more  directly  useful  for  stable  grasping. 

•  There  is  a  complex  relation  between  displacement  profile  and  surface  stress 
profile. (Phillips  and  Johnson,  1981) 

In  this  paper,  the  elastic  characteristics  of  the  material  are  assumed  to  be  ho¬ 
mogeneous  and  isotropic.  In  addition,  the  dimension  of  the  contact  area  will  be 
assumed  to  be  infinitely  large  in  a  direction,  say  s,  such  as  a  line  load.  By  that  as¬ 
sumption,  we  just  need  to  analyze  two  dimensional  behavior  in  a  slice  perpendicular 
to  that  direction. 

From  the  theory  of  elasticity (  Timoshenko  and  Goodier,  1951),  under  above 
assumptions  about  the  material,  the  differential  equation  of  the  equilibrium  and 
the  boundary  condition  for  the  stress  are  linear.  Consequently,  it  is  possible  to 
consider  the  interior  stress  or  strain  distribution  due  to  a  general  contact  as  the 
superposition  of  those  quantities  due  to  a  set  of  line  contacts. 

We  now  study  the  behavior  of  an  infinite  homogeneous  and  isotropic  elastic 
material  under  line  contact  with  negligible  contact  width.  For  a  two  dimensional 
problem  in  polar  coordinates,  the  distribution  of  stress  follows  a  simple  radial  dis¬ 
tribution.  From  the  theory  of  elasticity (  Timoshenko  and  Goodier,  1951),  for  the 
concentrated  force  inclined  from  the  vertical  by  an  angle,  a  <  jt/2,  we  have 
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2  Feo*(6  —  a) 

Or  ~  - 

or 
0$  —  0 

*rt  =  0 

where  <rr  is  the  radial  stress  at  the  point  (r,#),  ot  is  the  stress  in  the  plane  at  (r,0) 
norma]  to  the  radial  stress,  rrt  is  the  shearing  stress  in  the  r,#-plane,  F  is  the  force 
per  unit  length  and  r  is  the  distance  from  the  point  of  application.  Note  that  we 
have  defined  the  normal  stress  to  be  positive  when  it  produces  compression.  Fig.l. 
shows  the  above  variables. 


F 


Fig.  1.  Line  force  on  elastic  half»plane 
Applying  the  tensor  transformation 

om  *  o,eo»*$ 
ow  *  oT tin* 9 

we  get  the  stress  expression  in  Cartesian  coordinates 
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(1) 
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2  F 


om  =  — cos[a  —  9)003*0 


wr 


2F 


o9  =  —  cos  (a  —  0)sin20 


where  r  =  yjx1  +  y*,  eoe0  —  */r,  em0  =  y/r.  By  decomposing  F  into  vertical 
component,  Fv  ,  and  tangential  component,  Ft  ,  we  have 


*  I 

etna  =  - 

F„ 

cosa  =  — 

F 

Then  Eq.(l)  can  be  expressed  as 


(2) 


a-  =  ^T*y{xFv+vFt)  (3) 

°v==  *(z>  +  y*)^xF*  +  vFt) 

We  will  only  consider  ot  since,  in  y  direction,  it  converges  faster  than  ov  as 
y  — *  ioo.  Applying  the  principle  of  superposition,  we  have,  for  a  load  distributed 
along  the  y-axis. 


«.(*,»)  -  r  -  vo)f.(tto)  +  -  n,)F,(ito)]<(ito  (4) 

/-«o 


where 


2x* 

*(**  +  y*)* 
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v)  = 


2  x»y 

ff(xs  +  y*)2 


for  y  €  R  and  x  6  (0,  oo).  F„(y)  and  F<(y)  are  the  surface  stress  distributions  in 
x  and  y  direction,  respectively.  Fig.  2.  shows  the  JfJ‘rM,(x,y)  and  F/tr"*(x,y)  for 


x  =  1  and  x  =  2. 
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Fig.  2.  Functions  of  K*ri"  and  Kf™  at  x  —  1  and  x  =  2 
Since  many  pressure  sensors  have  outputs  based  on  the  strain  of  the  sensor  , 
which  is  the  fractional  change  in  the  linear  dimensions  of  a  small  cubic  volume 
element,  e.g.  electric-resistance  strain  gauges,  we  need  to  find  the  relations  between 
the  components  of  stresses  and  the  components  of  the  strains.  By  Hooke’s  Law  in 
the  theory  of  elasticity,  we  have 


««  =  o*  ~  "{o*  +  *.)] 

~  +  *«)]  (5) 
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where  E  is  the  modulus  of  elasticity,  v  is  Poisson’s  ratio,  e,  e„  and  e,  are  the 
components  of  strains  along  x,y,  z  direction  ,  respectively. 

Using  the  plane  strain  assumption  which  states  that  for  a  line  force  of  infinite 
extension  on  an  elastic  half-space,  the  strain  in  the  direction  of  the  line,  cs,  must 
be  tero  by  symmetry.  Then  Eq.(5)  becomes 

£* =  ~  •'V*  ~  (6) 

€v  =  -  v')°v  ~  + *)*»] 

From  Eq(6),  we  can  get  the  expression  of  oz  by  ez  and  f, 

Ot  =  Af(l  -  |/)C«  +  Uyty]  (7) 

where 

E 

A=  (1  +  i/)(l  —  2i/) 

It  should  be  note  that  if  */=0.5,  from  Eq.(6),  cE  =  -£,  =  ^(0,  —  o„).  That 
means  there  is  not  enough  information  from  (t  and  iv  to  determine  at  or  ay  .  In 
this  paper,  when  we  use  the  strain  distribution  in  x  direction  as  the  only  sensed 
data  to  solve  the  inverse  problem,  we  always  assume  that  v  is  0.5.  Actually,  the 
Poisson’s  ratio  of  rubberlike  materials,  which  usually  are  used  to  cover  robot  fingers 
and  tactile  sensors,  is  0.5  approximately. 

We  will  see  in  the  end  of  section  4,  the  surface  force  profile  reconstructed  from 
stress  distribution  is  more  accurate  than  from  strain  distribution  after  they  are 
truncated. 


Again,  using  the  principle  of  superposition  for  general  contact,  we  have 


«.(*,»)  =  /”  -  Vc)n(Vo)  +  Kf“'(x,V-  yo)F,(Vo)]</yo  (8) 

J—oo 

where 


zrtfrain / _  *  _  3x(x*  -  y1) 

K%  (  ’ V)  ~  2srJ5’(**  +  y*j* 

tfHrain/  \  _  3y(l  1/  ) 

K *  (I,V)  “  2V£(iH-^)> 


for  y  €  R  and  x  e  (0,oo).  Fig.  3.  shows  the  K[train{x,y)  and  Jcfj‘ro,n(x,y)  for  x  =  1 
and  x  =  2. 
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Fig.  3.  Functions  of  K[,rain  and  K[tr“4n  at  x  =  1  and  x  =  2. 

In  sake  of  convenience,  we  write  Eq.(4)  and  Eq.(8)  in  the  general  form  as 

g(y)  rn  f  (Jfi(y  -  yo)F.(yo)  +  Kt(y  -  yo)F,(yo)ldyo 

J  —00 


(®) 
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»\V  IV  • 
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where  y  e  R.  We  will  indicate  that  KitKj  and  g  are  in  the  form  of  stress  or  strain 
when  it  is  necessary.  In  Eq.(d),  there  is  no  x  dependence  since  we  assume  that  the 
array  of  tactile  sensors  is  usually  mounted  on  a  horizontal  plane  beneath  the  surface 
of  elastic  material.  Thus,  z  in  Eq.(9)  is  a  constant  which  expresses  the  depth  of  the 
sensors  from  the  surface. 
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3.  The  DFT  approach  for  inverse  problem 

In  this  section,  a  method  to  solve  Eq.(9)  approximately  will  be  examined.  We 
will  assume  that  Fv(y)  is  produced  by  the  strip  contact,  i.e.  it  satisfies 

Fv =  /  *i(»).  K  V  €  [«i,oj]; 

'  \0,  oterwise, 

where  Ji(y)  is  a  continuous  function  defined  on  [ai,a2]  ,  Fi(y)  >  0  for  y  €  (ai,a2) 

and  Fi(y)  >  0  at  y  =  a*  and  y  =  a2,  for  some  finite  a*  and  a2.  Moreover,  we  will 

assume  that  |  a,  |<  A/2,  i  =  1,2,  where  A  is  the  half  width  of  surface  of  the  finger. 

We  made  this  assumption  because,  even  for  a  very  narrow  force  profile,  a  broad 

distribution  of  stress  and  strain  will  be  produced  .  This  can  be  illustrated  by  Fig.2 

and  Pig.3  which  show  the  distribution  of  stress  and  strain  due  to  line  contact  when 

we  let  one  of  Fv  and  Ft  be  unit  impulse  function  and  another  be  zero.  In  other 

words,  for  a  strip  contact,  we  need  the  observations  which  distribute  more  broadly 

than  the  width  of  the  strip  to  obtain  the  surface  force  profile.  About  Ft(y ),  we  wi'l 

just  consider  two  special  cases  which  may  happen  during  the  grasping  process.  In 

addition,  we  will  assume  that  the  modulus  of  elasticity,  Et  is  one  in  the  rest  of  this 
paper. 

3.1  Surface  force  profile  without  tangential  component 

This  situation  may  happen  when  a  robot  hand  grasps  an  object  ,  which  lies  on 
a  stable  support,  to  determine  if  the  normal  force  is  large  enough  to  lift  the  object. 
In  this  case,  Eq.(9)  becomes 

g{y)=  f  Ki{y  -  yo)F,(yo)dyo  y€R  (10) 

OO 
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This  integral  equation  is  of  convolution  type,  in  which  g(y)  is  the  observation  ob¬ 
tained  by  sensors,  F„(y)  is  an  unknown  function,  Ky(y)  is  usually  called  the  kernel. 


Mathematically,  to  recover  the  unknown  function,  F„(y),  from  the  data  of  y(y) 
is  referred  to  as  the  deconvolution  problem  or,  in  general,  the  problem  of  solving  op¬ 
erator  equation  of  the  first  kind.  This  problem  is  ill-posed  in  the  sense  of  Hadamard 
(1923) .  Several  methods  for  this  problem  have  been  derived  which  all  try  to  find  a 
solution  by  supplementing  the  calculations  with  some  extra  information,  for  exam¬ 
ple,  generalized  inverse  and  regularization,  etc.(Hilgers,  1973). 

Next,  we  will  examine  a  method  depending  on  the  Fourier  Transform.  It  is 
known  that  the  Fourier  Transform  of  the  convolution  of  Ki  and  Fv  equals  to  the 
normal  multiplication  of  Fourier  Transform  of  K\  and  Fv  if 

/  I  *i(y)  I  dy  <  oo 

«/  — oo 

and 

/  I  -Mv)  \dy  <oo 

J— OO 

That  is, 


^HA(«)  =  iH  (11) 

where  Ki{u>),F„(u>)  and  g(u>)  are  Fourier  Transforms  of  Ki,Fv  and  g,  respectively. 

It  is  clear  that  if  K\{w)  ^  0,  V|u>|  <  oo  and  g[u)/Ki(u)  €  Lj(— oo.oo),  the  Eq.(10) 
can  be  solved  simply  as 


F.(V)  =  *"‘(^7)  V  6  * 

K  i(u>) 

11 


(12) 


Actually,  it  is  easy  to  check  that,  for  and  K[tTain  given  in  Eq.(4)  and 

Eq.(8),  their  Fourier  Transforms  are 

K[tr'“(ui)  =  e-H(l  +  x[u>\)  u€R 

and 


K?rain(u>)  =  1.5i|w|e",M  u  €  R 


Fig.  4.  The  Fourier  Transform  of  at  x  =  1.5. 


Fig.  5.  The  Fourier  Transform  of  K[tT*in  at  x  =  1.5  . 
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From  Fig.4.,  we  find  that  if  the  distribution  of  stress  can  be  obtained  directly 
or  indirectly  by  sensors,  the  inverse  problem  can  be  solved  by  Eq.(12).  If  we  have 
to  use  the  distribution  of  strain,  the  conditions  of  Eq.(12)  can  no  be  satisfied  since 
Kftra,n(0)  =  0,  (see  Fig. 5).  We  will  see  later  that  this  problem  will  be  solved  when 
the  observation,  g(y ),  is  truncated  and  DFT  is  applied. 

Physically,  the  continuous  expression  of  g[y)  can  not  be  obtained  directly,  but 
only  the  discrete  data  of  y(y).  So,  we  are  forced  to  apply  the  Discrete  Fourier 
Transform  to  solve  the  inverse  problem  approximately. 

The  discretization  of  Eq.  (10)  is 

siVn)  =  £  Kx(yn  -  V,)F,(y<) Ay,  n  =  0, ±1, ±2, . . .  (13) 

<«  — 00 

where  Ay,  =  yt+j  —  y,.  If  we  let  Ay,  be  identical  for  all  *,  denoted  by  A,  which 
will  be  called  the  sampling  period,  and  denote  y(y„)  by  y(n),  Ki(yn  —  y.) Ay,  by 
Ki (n  —  i)  and  F*(yi)  by  Fv(i),  then  Eq.(13)  can  be  written  as 

y(n)  =  22  n  =  0, ±1,±2, ...  (14) 

«s -00 

This  is  the  discrete  convolution  equation.  As  we  have  mentioned  before,  physi¬ 
cally,  it  is  impossible  to  get  the  data  of  observation,  y(n),  everywhere.  We  have  to 
cut  off  y(n),  at  n  —  -N  and  N, where  JVA  =  A.  Due  to  the  properties  of  K[tr,tt(y) 
and  K[lrain{y)  and  the  assumption  about  Fv(y),  if  A  is  large  enough,  g(N)  should 
be  very  small.  Then,  we  can  approximate  Eq.(14)  as 

g(n)  =  £  Ki(n  -  i)Fv{i)  n  =  0,±1,±2 . ±N  (15) 

i=-N 
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g(k)  =  Ki(k)Fv(k)  k  =  0,  ±1,  ±2, ....  ±N  (16) 

where  g(k)  =  9 (t)*~,J**k/aN+\  j  =  I,  *  =  0,±1,... ,±NAs  the  DFT  of 

y(i),  and  similar  to  Ki(k)  and  Fv(k).  Fig.6  and  Fig.7  give  the  DFT  of  K[trttt  and 

K[trmin.  I 

From  the  theory  of  digital  signal  processing  (Chen, 1079),  we  know  that  ,if  the 
Fourier  Transform  of  a  function  is  almost  zero  when  |u;|  is  greater  than  the  half  of  the 
sampling  frequency,  2x/A,  it  will  be  equal  to  the  discrete  Fourier  Transform  of  that 
function  times  sampling  period  at  u  =  2kir/((2N  +  1)A)  for  k  —  0,±1,±2,...,±N, 
where  A  is  the  sampling  period.  Comparing  Fig.4  and  Fig.5  with  Fig.6  and  Fig.7, 
we  find  K[tr,t,[y)  and  K[tT*in{y)  are  almost  satisfied  above  condition  and  result. 
But,  for  K[iraint  they  are  different  at  0.  The  reason  is  that  our  function  is  not 
periodic  and  N  can  not  be  arbitrary  large,  i.e.  it  have  to  be  truncated.  This  just 
give  us  a  chance  to  apply  the  DFT  to  solve  the  inverse  problem. 

Applying  Inverse  Discrete  Fourier  Transform,  the  unknown  function,  -F„(fc),  is 

F,(t)  “  wn  ‘-o-*1-**  <l7> 

In  next  section,  by  example,  we  will  find  that  the  reconstructed  function  pro¬ 
duced  by  Eq.(17)  is  very  accurate  with  respect  to  the  designed  one.  We  will  also 
discuss  the  effect  of  truncation  and  the  choice  of  K\. 

3.2  Surface  force  profile  with  constant  tangential  component 

This  situation  may  happen  when  a  robot  hand  lifts  an  object  which  has  a 
rectangular  surface.  In  order  to  avoid  slipping  during  grasping,  we  are  interested  in 
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the  tangential  force  on  the  average,  instead  of  at  every  point.  So,  we  assume  that 


Ft 


if  y  €  [ai,  a;]; 
otherwise, 


where  c  is  a  constant,  and  the  relation  of  Ft[y)  and  Fv(y)  at  the  boundary  of  cone 
of  friction  is  estimated  by 


[  Ft[y)dy  =  n  f  Fv(y)dy. 

J  —  oo  J  —00 

By  assumption  of  Ftt 

«=— —  rF"Wdv  (i8> 

aj  —  aj  /•, 

where  n  is  the  coefficient  of  static  friction.  Moreover,  we  assume  that  a  -  aj  =  a2 
since  we  will  need  to  determine  the  boundary  of  the  load  from  g[y).  Under  above 
assumptions,  Eq.(O)  becomes 

»(y)  =  J“jK ,(y  -  vo)F.(lfc)  +  £*,(»  - 1»)  /_”  F.(u)du}dy0  (19) 
Next,  we  will  discuss  how  to  solve  Eq.(l9)  numerically. 

By  discretizing  Eq.(l9)  and  truncating  the  observation,  we  have 


9<n)=  £  K.tn-OF.W  +  i^ff:  «(•)*.(»-))(  E  *(*)) 


im-N 


n  =  0,±1,...'±JV 

in  which  M  is  taken  such  that  2 a  =  2M A,  and  u(t)  is  defined  as 


(20) 


m_fl,  for  *:»A€  [-«,«); 
'  '  10,  otherwise. 
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Applying  DFT  to  Eq.(20),  we  have 


E  »(")“”  =  (  E  *.(»)«>“)(  E  F.(n)«.~) 

n=—N  n--N  ns-Af 


+  5T7(  £  ft(»))(  £  *>("Kr)(  £  tt(n)twnr) 

n=-N  n=-N  n=-N 

t  =  0,±1  ,...t±N 
where  w  =  e“,,*r/,N+1,j  =  y/-l.  When  r  =  0, 


E  »(»)  =  (  E  *>(»))<  E  *<»))  +  ,,(2w  ?)(  E  *»(»))(  £  *•.<»» 


nm—N  tl*—N 


na*-N 


n=-N 


Since  /Tj(y)  is  a  odd  function,  we  have 


£  *i(»)  =  o 

n=-W 


Therefore,  Eq.(2l)  becomes 


F=  E  ft(»)  =  ( E  »(«))( E  *.(«)>-* 

iw-W  tue—S  nm—N 


Discretizing  Eq.(18),  we  have 


Return  to  Eq.(20),  we  have 
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M  S 

g[n)  —  e  Ka(n  —  *)  =  K\(n  —  *)F„(t)  n  =  0,±1,...,±JV  (23) 

i=-M  x=-N 

Denoting  left  hand  side  of  above  equation  as  G(n),  Eq.(23)  becomes 

N 

G(n)  =  51  Ki(.n~  «W0  n  =  0,±1, ...  ±  N  (24) 

i=-N 

This  is  the  equation  we  have  solved  in  above  subsection.  However,  the  problem 
of  this  subsection  has  not  solved  since  we  do  not  know  M  yet.  As  we  have  assumed 
before,  F„(y)  and  Ft(y)  have  symmetric  support,  [-a,®].  And  we  will  see  in  next 
section  that  the  algorithm  in  subsection  3.1  has  &  very  good  accuracy.  By  these 
properties,  we  can  find  M .  The  method  is  shown  as  follows. 

We  have  known  that 

/OO 

Ki[y  -  yo)F„(y0)  +  Kt{y  -  yo)F,(y0)dy0 

-00 

Then, 

0{-V )  =  f  Ki[-V  -  VoWlfo)  +  Ki{-y  -  yo)Ft(yo)dyo 

»  —  OO 

By  changing  the  integral  variable  and  utilizing  the  odd  and  even  property  of 
Jfi(y)  and  Jf](y),  we  have 

f(-v)  =  /  Ki{v  ~  vo)^.(-vo)  -  *,(v  -  ife)fi(-vb)<*ve 

J-QO 

Since  we  assumed  that  Ft(y)  is  a  constant,  taking  the  average  of  g[y)  and  y(— y), 
we  have 
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y(v)  =  J  ^  Ki(y  -  y0)Fv{vo)dyo  (25) 

where  y(y)  =  (y(y)  +  y(-y)]/2,  F.(y)  =  [F,(y)  +  F*(-y)]/2. 

By  the  assumption  of  F*(y),  Fv(y)  still  have  a  symmetric  support,  [-a,  a].  There¬ 
fore,  applying  the  DFT  to  Eq.(25),  we  can  find  M  by  the  algorithm  in  subsection  3.1. 
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4.  Examples 

In  this  section  we  will  give  some  examples  for  some  special  load  profiles  and  com¬ 
pare  the  results  which  are  got  from  the  observations  expressed  as  the  distribution 
of  stress  and  distribution  of  strain  when  t/  =  0.5. 

Some  parameters  are  fixed  as  follows. 

•  the  depth  of  the  sensors,  x  =  1.5 

•  the  sampling  period,  A  =  0.5 

•  the  modulus  of  elasticity,  E  —  1 

•  the  numbers  of  elements  of  sensors,  2N  +  1  =  121 

We  will  consider  two  kinds  of  rigid  load:  cylinder  and  rectangle,  as  shown  in 
Fig.8. 


V 


|  elastic  material 
x 


V 

elastic  materia] 


x 


Fig.8  The  cylindrical  and  rectangular  indentor 
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From  Conway  et  al.  (1966),  for  cylinder  indenting  an  elastic  half-plane,  the 
surface  stress  is  given  by 


f  m  =  I  ~  for  y 6  [— °. a); 

vKy>  \0,  otherwise, 


where  a  is  the  half-width  of  the  contact  region  and  P  is  the  force  per  length;  for 
rectangular  indentor,  the  stress  on  the  surface  is  given  by 


*(v) 


=  (  V“a-va’ 

lo, 


for  y  6  [-a,e]; 
otherwise, 


For  numerical  simulation,  we  will  adjust  o  by  o  +  t  for  e  C  1 ,  since,  physically, 
F„(a)  can  not  be  infinite  large. 

In  following  examples,  g[y)  is  produced  by  the  discrete  form  of  Eq.(4)  and  Eq.(8) 
directly.  To  evaluate  the  reconstructed  surface  stress, Fr,  with  respect  to  designed 
one,  F,  we  use  the  absolute  error  defined  as 


e=uhi£K'F'(<)-m- 


(26) 


4.1  Cylinder  Indenting  without  friction 

The  distributions  of  stress  and  strain  at  x  =  1.5  produced  by  cylindrical  indentor 
are  given  in  Fig.9  and  Fig.ll,  respectively,  in  which  the  designed  surface  stresses 
are  also  displayed.  Fig.10  and  Fg.12  give  the  surface  stress  reconstructed  from  the 
observations  given  in  Fig.9  and  Fig.ll,  respectively.  For  comparing,  the  designed 
surface  stresses  are  displayed  again. 


By  Eq.(26),  the  errors  for  these  two  cases  are 

ettr‘“  =  3.93  x  10-6  ,e#‘r“n  =  4.04x1(Tb 

4.2  Rectangle  indenting  without  friction 

Fig.13-Fig.16  display  the  same  things  given  in  above  subsection  for  rectangular 
indenting. 


-1.0 


Fig.13  1.  the  designed  surface  stress;  2.  the  stress  at  z  =  1.5 


Ml  >1 


-1.0 


-3.0  -1.0  1.0  8.0 

Fig.14  1.  the  designed  surface  stress;  2.  the  reconstructed  surface  stress 
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-3.0  -1.0  1.0  3.0 

Fig.16  1.  the  designed  surface  stress;  2.  the  reconstructed  surface  stress 
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By  Eq.(26),  the  errors  for  these  two  cases  are 

e.lr.«  _  4A4  x  10-6.  j  Qg  x  jq-4 


4.3  Rectangular  indenting  with  constant  friction 

The  distributions  of  stress  and  strain  at  x  =  1.5  produced  by  cylindrical  indentor 
are  given  in  Fig.17  and  Fig. 19,  respectively,  in  which  the  designed  surface  stresses 
in  vertical  and  tangential  direction  are  also  displayed.  Fig. 18  and  Fig.20  give  the 
surface  stresses  reconstructed  from  the  observations  given  in  Fig.17  and  Fig. 19, 
respectively.  For  comparing,  the  designed  surface  stresses  are  displayed  again. 


Fig. 18  1.  the  designed  vertical  surface  stress; 

2.  the  reconstructed  vertical  surface  stress 

3.  the  designed  tangential  surface  stress 

4.  the  reconstructed  tangential  surface  stress 


-3.0  -1.0  1.0  3.0 

Fig.19  1.  the  designed  vertical  surface  stress; 

2.  the  designed  tangential  surface  stress;  3.  the  strain  at  x  =1.5 
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Fig. 20  1.  the  designed  vertical  surface  stress; 

2.  the  reconstructed  vertical  surface  Btress 

3.  the  designed  tangential  surface  stress 

4.  the  reconstructed  tangential  surface  stress 
By  Eq.(26),  the  errors  for  these  two  cases  are 

c**r‘M  =  6.21  x  10_s,  ej*""  =  1.15  x  10“* 

e,vtTain  =  2.95  x  10"*,  e?rmin  =  1.45  x  10~* 

where  the  subscripts  v  and  t  are  expressed  for  vertical  and  tangential  component, 
respectively. 

From  above  examples,  we  find  that  e*‘r*’n  >  e**r"*  for  all  cases.  There  are  many 
sources  to  produce  errors.  One  of  them  is  that  the  observation  has  been  truncated. 
We  have  known  that  the  further  the  point  is  from  the  area  of  surface  load,  the  little 
the  effect  of  surface  force  to  strain  or  stress  beneath  the  surface  is.  Comparing  the 
distributions  of  stress  and  strain  for  the  same  load,  e.g.  Fig.9  and  Fig.ll,  we  find 
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that  the  effect  to  the  strain  is  broader  than  the  stress.  Thus,  more  information 
will  be  lost  when  we  truncate  the  observation  which  is  expressed  as  strain  instead 
of  stress.  Above  argument  can  be  illustrated  by  Fig.21  and  Fig.22,  in  which  the 
number  of  element  of  sensors  is  only  41. 


imnwsii 


Fig. 21  the  surface  stress  reconstructed  from  the  distribution  of  stress  at  x 


-.100 


-1.0 


-.20 


Fig.22  the  surface  stress  reconstructed  from  the  distribution  of  strain  at  x  =  1.5 
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5.  Discussion 

In  this  paper,  we  described  how  to  express  the  distribution  of  stress  and  strain 
as  the  convolution  of  the  surface  stress  and  relative  kernel.  We  gave  the  DFT 
approach  for  two  special  cases.  By  example,  we  discussed  the  effect  of  truncation 
of  the  observation  and  the  effect  of  choosing  as  observation  the  strain  or  stress. 

As  we  have  mentioned  in  first  section,  there  are  several  methods  to  solve  op¬ 
erator  equations  of  first  kind.  By  comparing  the  DFT  approach  with  one  of  the 
regularization  methods  of  Hilgers(1973),  we  found  that  the  former  has  the  proper¬ 
ties  of  higher  computing  speed  and  better  accuracy.  However,  if  noise  with  high 
frequency  is  added  to  the  observation,  one  finds  that  the  DFT  method  is  much 
more  sensitive  to  the  noise,  even  when  the  deviation  of  the  noises  is  very  small. 
It  should  be  noted  that,  for  both  of  these  approaches,  the  analog  networks  have 
been  introduced  to  solve  the  inverse  problem  by  parallel  computation  (Poggio,1985; 
Peckerer  ).  This  should  be  of  interest  in  real  time  control  of  the  robot  fingers. 

There  are  several  problems  that  remain  to  be  solved.  One  is  how  to  reconstruct 
the  surface  stress  profile  when  it  has  a  general  tangential  distribution.  To  solve  this 
problem,  more  information  on  the  strain  or  the  stress  is  required.  Sensors  which 
are  mounted  on  two  horizontal  planes  beneath  the  surface  of  elastic  material  may 
be  helpful  for  this  case.  The  other  problem  is  how  to  deal  with  three  dimensional 
problem,  which  is  important  for  physical  situations.  To  solve  this  problem,  a  more 
complex  model  of  the  elastic  material  should  be  applied. 
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